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Abstract 

In the reaction e + e~ — > WW — > (qiq 2 )(q3q 4 ) the usual hadronization models 
treat the colour singlets qiq 2 and q^q^ coming from two W bosons indepen- 
dently. However, since the final state partons may coexist in space and time, 
cross-talk between the two evolving hadronic systems may be possible during 
fragmentation through soft gluon exchange. This effect is known as Colour 
Reconnection. In this article the results of the investigation of Colour Recon- 
nection effects in fully hadronic decays of W pairs in DELPHI at LEP are 
presented. Two complementary analyses were performed, studying the particle 
flow between jets and W mass estimators, with negligible correlation between 
them, and the results were combined and compared to models. In the frame- 
work of the SK-I model, the value for its k parameter most compatible with 
the data was found to be: 

«sk-i — 2.2+^3 

corresponding to the probability of reconnection V TCCO to be in the range 
0.31 < P r eco < 0.68 at 68% confidence level with its best value at 0.52. 
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1 Introduction 



The space-time development of a hadronic system is still poorly understood, and 
models are necessary to transform a partonic system, governed by perturbative QCD, to 
final state hadrons observed in the detectors. 

WW events produced in e + e~ collisions at LEP-2 constitute a unique laboratory to 
study and test the evolution of such hadronic systems, because of the clean environment 
and the well-defined initial energy in the process. Of particular interest is the possibility 
to study separately one single evolving hadronic system (one of the W bosons decaying 
semi-leptonically, the other decaying hadronically) , and compare it with two hadronic 
systems evolving at the same time (both W bosons decaying hadronically). 

Interconnection effects between the products of the hadronic decays of the two W 
bosons (in the same event) are expected since the lifetime of the W bosons (rw — ^/Tw — 
0.1 fm/c) is an order of magnitude smaller than the typical hadronization times. These 
effects can happen at two levels: 

• in the evolution of the parton shower, between partons from different hadronic sys- 
tems by exchanging coloured gluons [1] (this effect is called Colour Reconnection 
(CR) for historical reasons); 

• between the final state hadrons, due to quantum-mechanical interference, mainly 
due to Bose-Einstein Correlations (BEC) between identical bosons (e.g. pions with 
the same charge). 

A detailed study by DELPHI of this second effect was recently published [2]. 

The first effect, the possible presence of colour flow between the two W hadronization 
systems, is the topic studied in this paper. This effect is worthy of study in its own right 
and for the possible effects induced on the W mass measurement in fully hadronic events 
(see for instance [3] for an introduction and [4] for an experimental review). 

The effects at the perturbative level are expected to be small [3], whereas they may 
be large at the hadronization level (many soft gluons sharing the space-time) for which 
models have to be used to compare with the data. 

The most tested model is the Sjostrand-Khoze "Type 1" CR model SK-I [5]. This 
model of CR is based on the Lund string fragmentation phenomenology. The strings are 
considered as colour flux tubes with some volume, and reconnection occurs when these 
tubes overlap. The probability of reconnection in an event is parameterised by the value 
k, set globally by the user, according to the space-time volume overlap of the two strings, 



The parameter k was introduced in the SK-I model to allow a variation of the percentage 
of reconnected events and facilitate studies of sensitivity to the effect. In this model 
only one string reconnection per event was allowed. The authors of the model propose 
the value of k — 0.66 to give similar amounts of reconnection as other models of Colour 
Reconnection. By comparing the data with the model predictions evaluated at several 
k values, it is possible to determine the value of k most consistent with the data and 
extract the corresponding reconnection probability. Another model was proposed by the 
same authors, considering the colour flux tubes as infinitely thin, which allows for Colour 
Reconnection in the case the tubes cross each other and provided the total string length 
is reduced (SK-II'). This last model was not tested. 

Two further models are tested here, these are the models implemented in HERWIG [6] 
and ARIADNE [7] Monte Carlo programs. In HERWIG the partons are reconnected, with a 



overlap • 




(1) 
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reconnection probability of 1/9, if the reconnection results in a smaller total cluster mass. 
In ARIADNE, which implements an adapted version of the Gustafson-Hakkinen model [8], 
the model used [9] allows for reconnections between partons originating in the same W 
boson, or from different W bosons if they have an energy smaller than the width of the 
W boson (this model will be referred as 'AR-2'). 

Colour Reconnection has been previously investigated in DELPHI by comparing in- 
clusive distributions of charged particles, such as the charged-particle multiplicity dis- 
tribution or the production of identified (heavy) particles, in fully hadronic WW events 
and the distributions in semi-leptonic WW events. The investigations did not show any 
effect as they were limited by statistical and systematic errors and excluded only the 
most extreme models of CR (see [10]). 

This article presents the results of the investigations of Colour Reconnection effects in 
hadronically decaying W pairs using two techniques. The first, proposed by L3 in [11], 
looks at the particle flow between the jets in a 4-jet WW event. The second, proposed 
by DELPHI in [12], takes into account the different sensitivity to Colour Reconnection of 
several W mass estimators. The first technique is more independent of the model and it 
can provide comparisons based on data. The second technique is more dependent on the 
model tested, but has a much larger sensitivity to the models SK-I and HERWIG. Since 
the particle flow and W mass estimator methods were found to be largely uncorrelated a 
combination of the results of these two methods is provided. 

The paper is organised as follows. In the next section, the LEP operation and the 
components of the DELPHI detector relevant to the analyses are briefly described. In 
section 3 data and simulation samples are explained. Then both of the analysis methods 
discussed above are described and their results presented in sections 4 and 5. The com- 
bination of the results is given in section 6 and conclusions are drawn in the seventh and 
final section. 

2 LEP Operation and Detector Description 

At LEP-2, the second phase of the e + e~ collider at CERN, the accelerator was operated 
at centre-of-mass energies above the threshold for double W boson production from 1996 
to 2000. In this period, the DELPHI experiment collected about 12000 WW events 
corresponding to a total integrated luminosity of 661 pb _1 . About 46% of the WW 
events are WW — > qiq2q3q4 events (fully hadronic), and 44% are WW — > qiq^P, where 
t is a lepton (semi-leptonic). 

The detailed description of the DELPHI detector and its performance is provided 
in [13, 14]. A brief summary of the main characteristics of the detector important for the 
analyses follows. 

The tracking system of DELPHI consisted of a Time Projection Chamber (TPC), the 
main tracking device of DELPHI, and was complemented by a Vertex Detector (VD) 
closest to the beam pipe, the Inner and the Outer Detectors in the barrel region, and two 
Forward Chambers in the end caps. It was embedded in a 1.2 T magnetic field, aligned 
parallel to the beam axis. 

The electromagnetic calorimeter consisted of the High density Projection Chamber 
(HPC) in the barrel region, the Forward Electromagnetic Calorimeter (FEMC) and the 
Small angle Tile Calorimeter (STIC) in the forward regions, complemented by detectors to 
tag the passage of electron-positron pairs from photons converted in the regions between 
the FEMC and the HPC. The total depths of the calorimeters corresponded to about 18 
radiation lengths. The hadronic calorimeter was composed of instrumented iron with a 
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total depth along the shortest trajectory for a neutral particle of 6 interaction lengths, 
and covered 98% of the total solid angle. Embedded in the hadronic calorimeter were 
two planes of muon drift chambers to tag the passage of muons. The whole detector was 
surrounded by a further double plane of staggered muon drift chambers. 

For LEP-2, the DELPHI detector was upgraded as described in the following. 

Changes were made to some of the subdetectors, the trigger system [15], the run 
control and the algorithms used in the offline reconstruction of tracks, which improved 
the performance compared to the earlier LEP-1 period. 

The major changes were the extensions of the Vertex Detector (VD) and the Inner 
Detector (ID), and the inclusion of the Very Forward Tracker (VFT) [16], which increased 
the coverage of the silicon tracker to polar angles with respect to the z-axis 1 of 11° < 
9 < 169°. To further improve the track reconstruction efficiency in the forward regions 
of DELPHI, the tracking algorithms and the alignment and calibration procedures were 
optimised for LEP-2. 

Changes were also made to the electronics of the trigger and timing system which 
improved the stability of the running during data taking. The trigger conditions were 
optimised for LEP-2 running, to give high efficiency for 2- and 4-fermion processes in the 
Standard Model and also to give sensitivity to events which may have been signatures 
of new physics. In addition, improvements were made to the operation of the detector 
during the LEP operating states, to prepare the detector for data taking at the very start 
of stable collisions of the e + e~ beams, and to respond to adverse background from LEP 
when it arose. These changes led to an overall improvement in the efficiency for collecting 
the delivered luminosity from about 85% in 1995, before the start of LEP-2, to about 
95% at the end in 2000. 

During the operation of the DELPHI detector in 2000 one of the 12 sectors of the 
central tracking chamber, the TPC, failed. After 1 st September it was not possible to 
detect the tracks left by charged particles inside the broken sector. The data affected 
corresponds to around 1/4 of the data collected in 2000. Nevertheless, the redundancy of 
the tracking system of DELPHI meant that tracks passing through the sector could still 
be reconstructed from signals in any of the other tracking detectors. As a result, the track 
reconstruction efficiency was only slightly reduced in the region covered by the broken 
sector, but the track parameter resolutions were degraded compared with the data taken 
prior to the failure of this sector. 

3 Data and Simulation Samples 

The analyses presented here use the data collected by DELPHI in the years 1997 to 
2000, at centre-of-mass energies y/s between 183 and 209 GeV. The data collected in the 
year 2000 with the TPC working in full, with centre-of-mass energies from 200 to 208 
GeV and a integrated luminosity weighted average centre-of-mass energy of 206 GeV, were 
analysed together. Data acquired with the TPC with a broken sector, corresponding to a 
integrated luminosity weighted average centre-of-mass energy of 207 GeV, were analysed 
separately and included in the results presented here. 

The total integrated luminosity of the data sample is 660.8 pb _1 , and the integrated 
luminosity weighted average centre-of-mass energy of the data is 197.1 GeV. 

To compare with the expected results from processes in the Standard Model including 
or not including CR, Monte Carlo (MC) simulation was used to generate events and 

1 The DELPHI coordinate system is a right-handed system with the z-axis collincar with the incoming electron beam, 
and the x axis pointing to the centre of the LEP accelerator. 
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simulate the response of the DELPHI detector. These events were reconstructed and 
analysed with the same programs as used for the real data. 

The 4-fermion final states were generated with the code described in [17], based on 
WPHACT [18], for the WW signal (charged currents) and for the ZZ background (neutral 
currents), after which the events were fragmented with PYTHIA [19] tuned to DELPHI 
data [20]. The same WW events generated at 189, 200 and 206 GeV were also fragmented 
with PYTHIA implementing the SK-I model, with 100% reconnection probability. The 
systematic effects of fragmentation were studied using the above WW samples and WW 
samples generated with WPHACT and fragmented with either ARIADNE [7] or HERWIG [6] at 
183, 189, 200 and 206 GeV. For systematic studies of Bose-Einstein Correlations (BEC), 
WW samples generated with WPHACT and fragmented with PYTHIA implementing the BE 32 
model [21] of BEC, were used at all energies, except at 207 GeV. The integrated luminosity 
of the simulated samples was at least 10 times that of the data of the corresponding year, 
and the majority corresponded to 100 times that of the data. 

To test the consistency of the SK-I model and measure the k parameter, large WW 
samples were generated in an early stage of this work with EXCALIBUR [22] at 200 and 
206 GeV, keeping only the fully hadronic decays. These samples were then fragmented 
with PYTHIA. It was verified for smaller subsets that the results using these large samples 
and the samples generated later with WPHACT are compatible. 

The qq(7) background events were generated at all energies with KK2f [23] and frag- 
mented with PYTHIA. For systematic studies, similar KK2f samples fragmented with 
ARIADNE [7] were used at 183, 189, 200 and 206 GeV. 

These samples will be referred to as "DELPHI samples" . 

At 189 GeV, to compare with the other LEP experiments and with different CR mod- 
els, 6 samples generated with K0RALW [24] for the 4-fermion final states were also used. 
These samples 2 will be referred to as "Cetraro samples" . The events in the different sam- 
ples have the final state quarks generated with the same kinematics, and differ only in the 
parton shower evolution and fragmentation. Three samples were fragmented respectively 
with PYTHIA, ARIADNE and HERWIG (using the tuning of the ALEPH collaboration), with 
no CR implementation. Three other samples were fragmented in the same manner but 
now implementing several CR models: the SK-I model with 100% reconnection proba- 
bility, the AR-2 model, and the HERWIG implementation of CR with 1/9 of reconnected 
events, respectively. 

4 The Particle Flow Method 

The first of the two analyses presented in this paper is based on the so-called "particle 
flow method". The particle flow algorithm is based on the selection of special event 
topologies, in order to obtain well defined regions between any two jets originating from 
the same W (called the Inside-W region) or from different Ws (called the Between-W 
region). It is expected that Colour Reconnection decreases (increases) particle production 
in the Inside-W (Between-W) region. Hence, by studying the particle production in the 
inter-jet regions it is possible to measure the effects of Colour Reconnection. However, 
this method requires a selection of events with a suitable topology (see below) which has 
a low efficiency (£25%). 

2 produced by ALEPH after the LEP-W Physics Workshop in Cetraro, Italy, October 2001 
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4.1 Event and Particle Selection 

Events with both Ws decaying into qiq2 are characterised by high multiplicity, large 
visible energy, and the tendency of the particles to be grouped in 4 jets. The background 
is dominated by qq(7) events. 

Charged particles were required to have momentum p larger than 100 MeV/ c and below 
1.5 times the beam energy, a relative error on the momentum measurement Ap/p < 1, 
and polar angle 9 with respect to the beam axis between 20° and 160°. To remove tracks 
from secondary interactions, the distance of closest approach of the extrapolated track 
to the interaction point was required to be less than 4 cm in the plane perpendicular to 
the beam axis and less than 4/sin# cm along the beam axis, and the reconstructed track 
length was required to be larger than 30 cm. 

Clusters in the electromagnetic or hadronic calorimeters with energy larger than 0.5 
GeV and polar angle in the interval 10° < 6 < 170°, not associated to charged particles, 
were considered as neutral particles. 

The events were pre-selected by requiring at least 12 charged particles, with a sum 
of the modulus of the momentum transverse to the beam axis, of charged and neutral 
particles, above 20% of the centre-of-mass energy. These cuts reduced the contributions 
from gamma-gamma processes and beam-gas interactions to a negligible amount. The 
momentum distribution of the charged particles for the pre-selected events is shown in 
Figure 1 and compared to the expected distribution from the simulation. A good agree- 
ment between data and simulation is observed. 




Figure 1: Momentum distribution for charged particles (range 0-50 GeV/c (a) and 
0-5 GeV/c (b)). Points represent the data and the histograms represent the contributions 
from simulation for the different processes (signal (white) and background contributions). 



About half of the e + e~— >qq(7) events at high-energy are associated with an energetic 
photon emitted by one of the beam electrons or positrons (radiative return events), thus 
reducing the energy available in the hadronic system to the Z mass. To remove these 
radiative return events, the effective centre-of-mass energy \fs', computed as described 
in [25], was required to be above 110 GeV. It was verified that this cut does not affect the 
signal from W pairs, but reduces significantly the contribution from the qq(7) process. 
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In the WW fully hadronic decays four well separated energetic jets are expected which 
balance the momentum of the event and have a total energy near to the centre-of-mass 
energy The charged and neutral particles in the event were thus clustered using the 
DURHAM algorithm [26], for a separation value of y cnt = 0.005, and the events were 
kept if there were 4 and only 4 jets and a multiplicity (charged plus neutral) in each 
jet larger than 3. The combination of these two cuts removed most of the semi-leptonic 
WW decays and the 2-jet and 3-jet events of the qq(7) background. The charged-particle 
multiplicity distribution for the selected events at 189 GeV is given in Figure 2, with data 
points compared to the histogram from simulation of signal and background processes. 
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Figure 2: Uncorrected charged-particle multiplicity distribution at a centre-of-mass en- 
ergy of 189 GeV. Points represent the data and the histograms represent the contribution 
from simulation for the different processes. 



For the study of the charged-particle flow between jets, the initial quark configuration 
should be well reconstructed with a good quark-jet association. At 183 GeV and above, 
the produced W bosons are significantly boosted. This produces smaller angles in the 
laboratory frame of reference between the jets into which the W decays, when compared 
to these angles at threshold (back-to-back). Hence, this property tends to reduce the 
ambiguity in the definition of the Between-W and Inside-W regions. The selection criteria 
were designed in order to minimize the situation of one jet from one W boson appearing 
in the Inside-W region of the other W boson. 

The selection criteria are based on the event topology, with cuts in 4 of the 6 jet-jet 
angles. The smallest and the second smallest jet-jet angle should be below 100° and not 
adjacent (not have a common jet). Two other jet-jet angles should be between 100° and 
140° and not adjacent (large angles). 

In the case that there are two different combinations of jets satisfying the above criteria 
for the large angles, the combination with the highest sum of large angles is chosen. This 
selection increases the probability to have a correct pairing of jets to the same W boson. 
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Table 1: Centre-of-mass energy (^/s in GeV), integrated luminosity (£ in pb -1 ), efficiency 
and purity of the data samples, number of selected events, number of expected events 
from 4-jet WW and background processes (total and separated by process), and efficiency 
of correct pairing of jets to the same W boson. 

The integrated luminosity, the efficiency to select 4-jet WW events and the purity of 
the selected data samples, estimated using simulation, and the number of selected events 
are summarised for each centre-of-mass energy in Table 1. The numbers of expected 
events are also given separately for the signal and the background processes, and were 
estimated using simulation. The efficiency to select the correct pairing of jets to the 
same W boson, estimated with simulation as the fraction of WW events for which the 
selected jets 1 and 2 (see later) correspond indeed to the same W boson, is given in the 
last column of the Table. 

The efficiency of the event selection criteria decreases with increasing centre-of-mass 
energy. This is primarily due to the 'large' angles being reduced as a result of the 
increased boost (becoming lower than the cut value of 100°) and 'small' angles being 
increased due to the larger phase-space available (becoming higher than the cut value of 
100°). Much for the same reason, the efficiency to assign two jets to the same W boson in 
the selected events increases slightly with increasing centre-of-mass energy, in opposition 
to what would happen at threshold with the W boson decaying into two back-to-back 
jets, that would never be selected to come from the same W boson by the requirement 
that their interjet angle should be between 100° and 140°. 

In the following analysis the jets and planar regions are labeled as shown in Figure 3: 
the planar region corresponding to the smallest jet-jet angle is region B in the plane 
made by jets 2 and 3; the second smallest jet-jet angle corresponds to the planar region D 
between jets 1 and 4 in the plane made by these two jets; the planar region corresponding 
to the greatest of the large jet-jet angles in this combination is region A and spans 
the angle between jets 1 and 2 in the plane made by these jets; and finally region C 
corresponds to the planar region spanned by the second large angle, between jets 3 and 
4 in the plane made by these two jets. In general, the planar regions are not in the same 
plane, as the decay planes of the W bosons do not coincide, and the large angles in this 
combination are not necessarily the largest jet-jet angles in the event. 

The distribution of the reconstructed masses of the jet pairings (1,2) and (3,4), after 
applying a 4C kinematic fit requiring energy and momentum conservation, is shown in 
Figure 4 (two entries per event). In the figure, data at 189 GeV (points) are compared 
to the expected distribution from the 4-jet WW signal without CR, plus background 
processes, estimated using the simulation (histograms). The contribution from the 4-jet 
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Figure 3: Schematic drawing of the angular selection. 



WW signal simulation is split between the case in which the two pairs of jets making the 
large angles actually come from their parent W bosons and the case in which the jets of 
a pair come from different W bosons (mismatch). 

4.2 Particle Flow Distribution 

The particle flow analysis uses the number of particles in the Inside-W and the 
Between-W regions. An angular ordering of the jets is performed as in Figure 3. The 
two large jet-jet angles in the event are used to define the Inside-W regions, and the two 
smallest angles span the Between-W regions, the regions between the different Ws. 

In general, the two W bosons will not decay in the same plane, and this must be 
accounted for when comparing the particle production in the Inside-W and Between-W 
regions. So, for each region (A, B, C and D) the particle momenta of all charged particles 
are projected onto the plane spanned by the jets of that region: jets 1 and 2 for region 
A; jets 2 and 3 for region B; jets 3 and 4 for region C; jets 4 and 1 for region D. Then, 
for each particle the rescaled angle $ reS caied is determined as a ratio of two angles: 



when the particle momentum is projected onto the plane of the region r. The angle $j is 
then the angle between the projected particle momentum and the first mentioned jet in 
the definition of the regions given above. The angle $ r is the full opening angle between 
the jets. Hence $ rescaled varies between and 1 for the particles whose momenta are 
projected between the pair of jets defining the plane. 

However, due to the aplanarity of the event about 9% of the particles in the data and in 
the 4-jet WW simulation have projected angles outside all four regions. These particles 
were discarded from further analysis. In the case where a particle could be projected 
onto more than one region, with < $ reS caied < 1, the solution with the lower momentum 
transverse to the region was used. This happened for about 13% of the particles in data, 
after background subtraction, and in the 4-jet WW simulation. 

This leads to the normalised particle flow distribution shown in Figure 5 at 189 GeV, 
where the rescaled angle of region A is plotted from to 1, region B from 1 to 2, region 



^rescaled ^i/^ 



(2) 
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Figure 4: Reconstructed dijet masses (after a 4C kinematic fit) for the selected pairs at 
189 GeV (2 entries per event) (see text). 

C from 2 to 3 and region D from 3 to 4. The statistical error on the bin contents (the 
average multiplicity per bin of $ re scaied divided by the bin width) was estimated using 
the Jackknife method [27], to correctly account for correlations between different bins. 
In this distribution the regions between the jets coming from the same W bosons (A and 
C), and from different W bosons (B and D), have the same scale and thus can be easily 
compared. 

After subtracting bin-by-bin the expected background from the observed distributions, 
we define the Inside-W (Between-W) particle flow as the bin-by-bin sum of regions A and 
C (B and D). These distributions are compared by performing the bin-by-bin ratio of 
the Inside-W particle flow to the Between-W particle flow. This ratio of distributions is 
shown for 189 GeV and 206 GeV in Figure 6. The data points are compared to several 
fully simulated WW MC samples with and without CR. 

A good agreement was found between the predictions using the WPHACT WW MC sam- 
ples and the predictions based on the KORALW WW MC samples, both for the scenario 
without CR and for the scenario with CR (SK-I model with 100% probability of recon- 
nection). For both sets of predictions the regions of greatest difference between the two 
scenarios span the rescaled variable $ reS caied from 0.2 to 0.8. 

4.3 Particle Flow Ratio 

After summing the particle flow distributions for regions A and C, and regions B and 
D, the resulting distributions are integrated from 0.2 to 0.8. The ratio R of the Inside-W 
to the Between-W particle flow is then defined as (with $ being the rescaled variable 

^rescaled) • 
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Figure 5: Normalised charged-particle flow at 189 GeV. The lines correspond to the sum 
of the simulated 4-jet WW signal with the background contributions (estimated from 
DELPHI MC samples), normalised to the total number of expected events (N events ). The 
dashed histogram corresponds to the sum with the simulated 4-jet WW signal generated 
by WPHACT with 100% SK-I. 
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Figure 6: The ratio of the particle flow distributions (A+C) / (B+D) at 189 GeV (a) and at 
206 GeV (b). The data (dots) are compared to WW MC samples generated with WPHACT 
(DELPHI samples) and K0RALW (Cetraro samples), both without CR and implementing 
the SK-I model with 100% probability of reconnection. The lines corresponding to WPHACT 
are hardly distinguishable from the lines corresponding to K0RALW in the same condition 
of implementation of CR. 
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Table 2: Values of the ratio R for each energy (errors are statistical only), and expected 
values with errors due to limited statistics of the simulation, all from DELPHI WPHACT 
WW samples. 



MC Sample 


X 2 /DF 


a, A 


[3,B 


7 


no CR 


7.31/5 


1.001 ±0.003 


(3.20 ± 0.36) x 10~ 3 


(-1.35 ±0.40) x 10" 4 


SK-I 100% 


1.46/1 


0.880 ±0.003 


(1.68 ± 0.44) x 10~ a 





Table 3: Results of the fit to the evolution of R with (^(GeV) — 197.5). 



R f™dn ch /d<S>(A + C)d$ 
f™dn ch /d$(B + D)d<S> ' 

To take into account possible statistical correlations between particles in the Inside-W 
and Between-W regions, the statistical error on this ratio R was again estimated through 
the Jackknife method [27]. 

The values for R obtained for the different centre-of-mass energies are shown in Table 2, 
and compared to the expectations from the DELPHI WPHACT WW samples without CR 
and implementing the SK-I model with 100% reconnection probability. These values for 
data and MC are plotted as function of the centre-of-mass energy in Figure 7. 

The changes in the value of R for the MC samples are mainly due to the dif- 
ferent values of the boost of the W systems. In order to quantify this effect 
a linear function R(^ - 197.5) = A + B • (y/s - 197.5) was fitted to the MC points 
with CR (with \/s in GeV), while for the points without CR the quadratic func- 
tion R(^/s — 197.5) = a + (3 ■ (y/s — 197.5) ± 7 • (y/s — 197.5) 2 was assumed (with y/s in 
GeV), giving reasonable x 2 /d.o.f. values. The fits yielded the results shown in Table 3. 

The MC without CR shows a stronger dependence on y/s. The function fitted to this 
sample was used to rescale the measured values of R for the data collected at different 
energies to the energy of 189 GeV, the centre-of-mass energy at which the combination 
of the results of the LEP experiments was proposed in [4]. All the rescaled values were 
combined with a statistical error-weighted average. The average of the R ratios rescaled 
to 189 GeV was found to be 



(R) = 0.979 ± 0.032(stat). 



(4) 
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Figure 7: The ratio R as function of y/s for data and MC (DELPHI WPHACT WW samples) , 
and fits to the MC with and without CR, and the combined ratio after rescaling all values 
to \fs = 189 GeV (see text). The value of the combined ratio at 189 GeV is shown at 
a displaced energy (upwards by 1 GeV) for better visibility, as well as all the values 
for the MC WW no CR' points and the corresponding fitted curve which are shown at 
centre-of-mass energies shifted downwards by 0.5 GeV. All errors for the MC values are 
smaller than the size of the markers. 

Performing the same weighted average when using for the rescaling the fit to the MC 
with CR, one obtains: 

(^CRrescaic) = 0.987 ± 0.032(stat). (5) 
Repeating the procedure, but now without rescaling the R ratios, the result is: 

(Korescale) = 0.999 ± 0.033(stat). (6) 

4.4 Study of the Systematic Errors in the Particle Flow 

The following effects were studied as sources of systematic uncertainties in this anal- 
ysis. 



4.4.1 Fragmentation and Detector response 

A direct comparison between the particle flow ratios measured in fully hadronic data 
and MC samples, i?4 q Data and -R4 q MC, respectively, is hampered by the uncertainties as- 
sociated with the modelling of the WW fragmentation and the detector response. These 
systematic uncertainties were estimated using mixed semi-leptonic events. In this tech- 



13 



nique, two hadronically decaying W bosons from semi-leptonic events were mixed together 
to emulate a fully hadronic WW decay. 

Mixing Technique 

Semi-leptonic WW decays were selected from the data collected by DELPHI at centre- 
of-mass energies between 189 and 206 GeV, by requiring two hadronic jets, a well isolated 
identified muon or electron or, in case of a tau candidate, a well isolated particle, all 
associated with missing momentum (corresponding to the neutrino) pointing away from 
the beam pipe. A neural network selection, developed in [28], was used to select the 
events. The same procedure was applied to the WPHACT samples fragmented with PYTHIA 
and HERWIG at centre-of-mass energies of 189, 200 and 206 GeV and with ARIADNE at 189 
and 206 GeV. The background to this selection was found to be of negligible importance 
in this analysis. Samples of mixed semi-leptonic events were built separately at each 
centre-of-mass energy for data and Monte Carlo semi-leptonic samples, following the 
mixing procedure developed in [2]. 

In each semi-leptonic event, the lepton (or tau-decay jet) was stripped off and the 
remaining particles constituted the hadronically decaying W boson. Two hadronically 
decaying W bosons were then mixed together to emulate a fully hadronic WW decay. 
The hadronic parts of W bosons were mixed in such a way as to have the parent W 
bosons back-to-back in the emulated fully hadronic WW decay. To increase the statistics 
of emulated events, and profiting from the cylindrical symmetry of the detector along 
the z axis, the hadronic parts of W bosons were rotated around the z axis, but were not 
moved from barrel to forward regions or vice-versa, as detailed in the following. 

When mixing the hadronic parts of different W events it was required that the two Ws 
had reconstructed polar angles back-to-back or equal within 10 degrees. In the latter case, 
when both Ws are on the same side of the detector, the z component of the momentum 
is sign flipped for all the particles in one of the Ws. 

The particles of one W event were then rotated around the beam axis, in order to have 
the two Ws also back-to-back in the transverse plane. Each semi-leptonic event was used 
in the mixing procedure between 4 and 9 times, to minimize the statistical error on the 
particle flow ratio R measured in the mixed semi-leptonic data sample. 

The mixed events were then subjected to the same event selection and particle flow 
analysis used for the fully hadronic events. The particle flow ratios -R m ixcd sl Data and 
-RmixedSLMC were measured in the mixed semi-leptonic data and MC samples, respectively, 
and are plotted as function of the centre-of-mass energy in Figure 8. 

The values of -R m i xe dSL measured in MC show a dependence on y/s. This effect is 
quantified by performing linear fits to the points measured with PYTHIA, ARIADNE and 
HERWIG, respectively. The differences between the measured slopes were found to be small. 
The function fitted to the PYTHIA points was used to rescale the values of R measured in 
data at different energies to 189 GeV. The rescaled values were then combined using as 
weights the scaled statistical errors. The weighted average R at 189 GeV for the mixed 
semi-leptonic events built from data was found to be 

( -Rmixcd SL Data ) 

1.052 ±0.027(stat). (7) 

For each MC sample, the ratio -Rmixcd SLData/-R m ixcdSLMC wa s used to calibrate the 
particle flow ratio measured in the corresponding fully hadronic sample, -R4 q MC; to 
compare it to the ratio measured in the data, (it^qData)- The correction factor 
-RmixedSL Data/-RmixcdSLMC was computed from the values of R rescaled to 189 GeV, calcu- 
lated from the fits to the mixed semi-leptonic samples built from the data and the MC. 
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Figure 8: The ratio -RmixcdSL as function of \fs for data and MC, and fits to the MC (see 
text). The ARIADNE points at 189 GeV and at 206 GeV have their centre-of-mass energy 
shifted and the error bars on data are tilted for readability. 



MC sample 


PYTHIA 


ARIADNE 


HERWIG 


R mixed SL Data / R mixed SL MC 
pCalibrated 
n 4qMC 


1.053 
1.018 


1.044 
1.011 


0.997 
1.004 



Table 4: Ratio of data to MC fitted values of R in mixed semi-leptonic samples, used 
to calibrate the i?4 q MC values for different models (upper line), and calibrated values of 
-R4qMO All values were computed at y/s = 189 GeV. 

The values for -RmixedSLMC are presented in Table 4, for the different models, along with 
the calibrated values of R^ q for the same models. 

The calibration factors differ from unity by less than 6%, and the largest difference of 
the calibrated i?4 q MC values when changing the fragmentation model, 0.014, was consid- 
ered as an estimate of the systematic error due to simulation of the fragmentation and of 
the detector response, and was added in quadrature to the systematic error. The error 
in the calibrated i?4 q MC values due to the statistical error on (-R m i xe dSLData) value used 
for the calibration, 0.026, was also added in quadrature to the systematic error. 

4.4.2 Bose-Einstein Correlations 

Bose-Einstein correlations (BEC) between identical pions and kaons are known to 
exist and were established and studied in Z hadronic decays in [29]. They are expected 
to exist with a similar behaviour in the W hadronic decays, and this is studied in [2]. 
They are implemented in the MC simulation samples with BEC via the BE 32 model 
of LUBOEI [21], which was tuned to describe the DELPHI data in [2]. However, the 
situation for the WW (ZZ) fully hadronic decays is not so clear, i.e. whether there are 
correlations only between pions and kaons coming from the same W(Z) boson or also 
between pions and kaons from different W(Z) bosons. The analyses of Bose-Einstein 
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correlations between identical particles coming from the decay of different W bosons do 
not show a significant effect [30] for three of the LEP experiments, whereas for DELPHI, 
an effect was found at the level of 2.4 standard deviations [2]. Thus, a comparison was 
made between the WPHACT samples without CR and with BEC only between the identical 
pions coming from the same W boson (BEC only inside), to the samples without CR 
and with BEC allowed for all the particles stemming from both W bosons, implemented 
with the BE 32 variant of the LUBOEI model (BEC all). The R values were obtained 
at each centre-of-mass energy, after which a linear fit was performed for each model to 
obtain a best prediction at f 89 GeV. The fit values were found to be in agreement to the 
estimate at 189 GeV alone, and for simplicity this estimate was used. The measurement 
of BEC from DELPHI of 2.4 standard deviations above zero (corresponding to BEC only 
inside), was used to interpolate the range of 4.1 standard deviations of separation between 
BEC only inside and BEC all. To include the error on the measured BEC effect, one 
standard deviation was added to the effect before the interpolation. The difference in the 
estimated values of R at y/s = 189 GeV, between the model with BEC only inside and 
the model with partial BEC all (at the interpolated point of 3.4/4.1), -0.013, was added 
in quadrature to the systematic error. 

4.4.3 qq(7) Background Shape 

The fragmentation effects, in the shape of the qq(7) background, were estimated by 
comparing the values of R obtained when the subtracted qq(7) sample was fragmented 
with ARIADNE instead of PYTHIA at the centre-of-mass energy of 189 GeV, and the differ- 
ence, 0.003, was added in quadrature to the systematic error. 

4.4.4 qq(7) and ZZ Background Contribution 

At the centre-of-mass energy of 189 GeV, the qq(7) cross-section in the 4-jet region 
is poorly known, due to the difficulty in isolating the qq(7) — > 4-jet signal from other 
4-jet processes such as WW and ZZ. The study performed in [31] has shown that the 
maximal difference in the estimated qq(7) background rate is 10% coming from changing 
from PYTHIA to HERWIG as the hadronization model, with the ARIADNE model giving 
intermediate results. Conservatively, at each centre-of-mass energy a variation of 10% 
on the qq(7) cross-section was assumed, and the largest shift in R, 0.011, was added in 
quadrature to the systematic error. 

The other background process considered is the Z pair production. The Standard 
Model predicted cross-sections are in agreement with the data at an error level of 10% [32]. 
The cross-section was thus varied by ±10% at each energy and the effect in R was found 
to be negligible. 

4.4.5 Evolution of R with Energy 

The R ratios were rescaled to y/s = 189 GeV using the fit to the MC without CR, 
however the correct behaviour might be given by the MC with CR. Hence, the difference of 
0.009 between the R values obtained using the two rescaling methods, using MC without 
CR (R) and with CR (.Rcr rescaie); w ^s added in quadrature to the systematic error. 

4.5 Results of the Particle Flow Analysis 

The final result for the average of the ratios R rescaled to 189 GeV is 
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PYTHIA no CR 


1.037 ±0.004 


PYTHIA SK-I 100% 


0.917 ±0.003 


ARIADNE no CR 


1.053 ±0.004 


ARIADNE AR2 


1.021 ±0.004 


HERWIG no CR 


1.059 ±0.004 


HERWIG 1/9 CR 


1.040 ±0.003 



Table 5: R ratios for the Cetraro samples at 189 GeV, calibrated with the mixed semi-lep- 
tonic events. 



(R) = 0.979 ± 0.032(stat) ± 0.035(syst). (8) 

In order to facilitate comparisons between the four LEP experiments, this value can 
be normalised by the one determined from simulation samples produced with the full 
detector simulation and analysed with the same method. The LEP experiments agreed 
to use for this purpose the Cetraro PYTHIA samples. These events were generated with the 
ALEPH fragmentation tuning but have been reconstructed with the DELPHI detector 
simulation and analysed with this analysis. The values of the R ratios obtained from the 
Cetraro samples at 189 GeV, calibrated using the mixed semi-leptonic events from these 
samples, are given in Table 5. 

The value of (R) measured from data is between the expected R ratios from PYTHIA 
without CR and with the SK-I model with 100% fraction of reconnection. The error 
of this measurement is larger than the difference between the values of R from ARIADNE 
samples without and with CR, and than the difference between values of R from the 
HERWIG samples without CR and with 1/9 of reconnected events. 

The following normalised ratios are obtained for the sample without CR and imple- 
menting the SK-I model with 100% CR probability, respectively: 



r toCR = p = °- 944 ± 0.031(stat) ± 0.034(syst), (9) 

-Rno CR 

r data = CR^data = ± Q (335^) ± .039(syst). (10) 

RcR 

In the above expressions, the statistical errors in the MC predicted values were propagated 
and added quadratically to the systematic errors on the ratios. 

It is also possible to define the following quantity, taking the predictions for Rcr and 
-RnoCR at y/s = 189 GeV from the PYTHIA samples in Table 5, 

Sr = ( ^ data) ~ p RnoCR = 0.49 ± 0.27(stat) ± 0.29(syst) , (11) 
Rcr — Rno cr 

from which it can be concluded that the measured (-Rdata) is compatible with intermediate 
probability of CR, and differs from the CR in the SK-I model at 100% at the level of 
1.3 standard deviations. The ability to distinguish between these two models can be 
computed from the inverse of the sum in quadrature of the statistical and systematic 
errors; it amounts to be 2.5 standard deviations. In Figure 9 the result of Sr is compared 
to the predicted values, in the scope of the SK-I model, as a funtion of the fraction of 
reconnected events. 
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Figure 9: Comparison of the measurement of the 5r observable to the predictions from 
the SK-I model as a function of the fraction of reconnected events. 

The result for the value of (R) can also be used to test for consistency with the SK-I 
model as a function of k and a log-likelihood curve was obtained. This also facilitates 
combination with the result obtained in the analysis in the following section, and for this 
reason the value of (R) is rescaled with PYTHIA without CR to a centre-of-mass energy 
of 200 GeV: the value obtained at 200 GeV is (i2)(200 GeV) = 1.024 ± 0.050. The 
values obtained for the predicted ratios Rn at 200 GeV and the log-likelihood curve, as 
a function of k, are shown in Figure 10. The value of k most compatible with the data 
within one standard deviation is 

* SK -t=4.13^ 7 . (12) 



5 Different Mw Estimators as Observables 

It has been shown [12] that the Mw measurement inferred from hadronically decaying 
W + W~ events at LEP-2, by the method of direct reconstruction, is influenced by CR 
effects, most visible when changing the value of k in the SK-I model. For the Mw(4q) 
estimator within DELPHI this is shown in [33]. Other published Mw estimators in LEP 
experiments are equally sensitive to k [34]. 

To probe this sensitivity to CR effects, alternative estimators for the Mw measure- 
ment were designed which have different sensitivity to k. In the following, the standard 
estimator and two alternative estimators, studied in this paper, are presented. The stan- 
dard estimator corresponds to that previously used in the measurement of the W mass 
by DELPHI [33]. Note that in the final DELPHI W mass analysis [35] results are given 
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Figure 10: a) Estimated ratio Rn at 200 GeV plotted as a function of different k values 
(top scale), or as function of the corresponding reconnection probabilities (bottom scale), 
compared to (R) measured from data after rescaling to 200 GeV (horizontal lines marked 
with R for the value and with la (2a) for the (R) value added/subtracted by one(two) 
standard deviations); the last three marks on the x axis, close to 100% of reconnection 
probability, correspond respectively to the values k = 100, 300, 800; b) corresponding 
log-likelihood curve for the comparison of the estimated values (Rn) with the data ((R)). 
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for the standard and hybrid cone estimators, with the hybrid cone estimator used to 
provide the primary result. The data samples, efficiencies and purities for the analysis 
corresponding to the standard estimator are provided in [33,35]. 

• The standard Mw estimator : 

This estimator is described in [33] and was optimised to obtain the smallest sta- 
tistical uncertainty for the W mass measurement. It results in an event-by-event 
likelihood Lj(Mw) for the parameter Mw- 

• The momentum cut Mw estimator : 

For this alternative Mw estimator the event selection was performed in exactly 
the same way as for the standard Mw estimator. The particle-jet association was 
also taken from this analysis. However, when reconstructing the event for the Mw 
extraction a tighter track selection was applied. The momentum and energy of the 
jets were calculated only from those tracks having a momentum higher than a certain 
Pcut value. An event-by-event likelihood L? cut (Mw) was then calculated. 

• The hybrid cone Mw estimator : 

In this second alternative Mw estimator the reconstruction of the event is the same 
as for the standard analysis, except when calculating the jet momenta used for the 
M w extraction. 




Figure 11: Illustration of the iterative cone algorithm within a predefined jet as explained 
in the text. 



An iterative procedure was used within each jet (defined by the clustering algorithm 
used in the standard analysis) to find a stable direction of a cone excluding some 
particles in the calculation of the jet momentum, illustrated in Figure 11. Starting 
with the direction of the original jet p^ d , the jet direction was recalculated (direction 
(1) on the Figure) only from those particles which have an opening angle smaller 
than R cone with this original jet. This process was iterated by constructing a second 
cone (of the same opening angle) around this new jet direction and the jet direction 
was recalculated again. The iteration was continued until a stable jet direction plf ne 
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was found. The jet momenta obtained, Pco ne , were rescaled to compensate for the 
lost energy of particles outside the stable cone, 

met 

f?J et _^ nJ ct . (i 1 *) 

fcone -Fcone rpjet ' V / 

cone 

The energies of the jets were taken to be the same as those obtained with the 
standard clustering algorithm (E jet conc — > E jct ). This was done to increase the 
correlation of this estimator with the standard one. The rescaling was not done 
for the p cut estimator as it will be used in a cross-check observable with different 
systematic properties. Again the result is an event-by-event likelihood Lp cono (M\y). 

Each of these previously defined Mw likelihoods had to be calibrated. The slope of 
the linear calibration curve for the Mw estimators is tuned to be unity, therefore only 
a bias correction induced by the reconstruction method has to be applied. This bias 
is estimated with the nominal WPHACT Monte Carlo events and the dependence on the 
value of re is estimated with the EXCALIBUR simulation. It was verified for smaller subsets 
that the results using these large EXCALIBUR samples and the samples generated with 
WPHACT are compatible. Neglecting the possible existence of Colour Reconnection (CR) in 
the Monte Carlo simulation results in event likelihoods L;(M W | event without CR), while 
Lj(M\v (event with CR) are the event likelihoods obtained when assuming the hypothesis 
that events do reconnect (100% CR in the scope of the SK-I model). To construct the 
event likelihoods for intermediate CR (values of re larger than 0) the following weighting 
formula is used : 

L;(M w |re) = [l-7 ? i(K)]-L i (M w |event without CR) + 7^ (re) • Lj (M w | event with CR) (14) 

where 7^ (re) is defined in Equation 1. The combined likelihood is produced for 
the event sample; the calibrated values for M w (re) were obtained for different val- 
ues of re using the maximum likelihood principle. In Figure 12 the difference 
dM w (re) = M w (re) — M w (re = 0) or the influence of re on the bias of the M w estimator is 
presented as function of re. 

The uncertainty on this difference is estimated with the Jackknife method [27] to 
take the correlation between M w (re) and M w (re = 0) into account. It was observed from 
simulations that the estimators dependency on re, for re below about 5, was not signifi- 
cantly different in the centre-of-mass range between 189 and 207 GeV. Therefore in the 
determination of re the dependency at 200 GeV was taken as default for all centre-of- 
mass energies. This value of centre-of-mass energy is close to the integrated luminosity 
weighted centre-of-mass energy of the complete data sample, which is 197.1 GeV. 

When neglecting the information content of low momentum particles or when using 
the hybrid cone algorithm, the influence of Colour Reconnection on the Mw estimator is 
decreased. The dependence of the estimator to re is decreased when increasing the 
value of p cut or when working with smaller cone opening angles R con e- 



5.1 The Measurement of k 

The observed difference AMw(std,i) = Mw std — Mw* in the event sample, where i 
is a certain alternative analysis, provides a measurement of re. When both estimators 
M w std and M w l are calibrated in the same hypothesis of re, the expectation values of 
AM w (std, i) will be invariant under a change of p cut or R CO nc- 

When neglecting part of the information content of the events in these alternative Mw 
analyses, by increasing p cut or decreasing R con c, the statistical uncertainty on the value of 
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Figure 12: The difference dM w (/t) = M w (/c) — M w (k = 0) is presented as a function of 
k, for different estimators. The curve for the standard estimator is the curve 
at the top. The curves obtained with the hybrid cone analysis for different values of the 
cone opening angle, starting from the top with 1.00 rad down to 0.75 rad, 0.50 rad and 
0.25 rad are indicated with dotted lines. The curves obtained with the momentum cut 
analysis for different values of p cut , starting from the top with 1 GeV/c, down to 2 GeV/c 
and 3 GeV/c are dashed. The vertical line indicates the value of k preferred by the SK-I 
authors [5] and commonly used to estimate systematic uncertainties on measurements 
using e + e~ — > W + W~ — > qiq 2 q3q 4 events. 
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the M w estimator is increased. Therefore a balance must be found between the statistical 
precision on AM w (std, i) and the dependence of this difference to k in order to obtain 
the largest sensitivity for a k measurement. This optimum was found using the Monte 
Carlo simulated events and assuming that the data follow the k — hypothesis, resulting 
in the smallest expected uncertainty on the estimation of k. 

For the p cut analysis an optimal sensitivity was found when using the difference 
AM w (std, p cut ) with p cut equal to 2 GeV/c or 3 GeV/c. Even more information about 
k could be extracted from the data, when using the difference AM w (std, R conc ), which 
was found to have an optimal sensitivity around R con e — 0.5 rad. No significant im- 
provement in the sensitivity was found when combining the information from these two 
observables. Therefore the best measure of k using this method is extracted from the 
AM w (std, Rconc = 0.5 rad) observable. Nevertheless, the AM w (std,p cut = 2GeV/c) ob- 
servable was studied as a cross-check. 

5.2 Study of the Systematic Errors in the AM W Method 

The estimation of systematic uncertainties on the observables AMw(std, i) follows 
similar methods to those used within the Mw analysis. Here the double difference is a 
measure of the systematic uncertainty between Monte Carlo simulation ('MC') and real 
data ('DA'): 

A syst (MC, DA) = | [M w std (MC) - M w std (DA)] - [Mw^MC) - Mw^DA)] | (15) 

where % is one of the alternative Mw estimators. The systematic error components are 
described below and summarised in Table 6. 

5.2.1 Jet Reconstruction systematics with MLBZs 

A novel technique was proposed in [36] to study systematic uncertainties on jet recon- 
struction and fragmentation in W physics measurements with high statistical precision 
through the use of Mixed Lorentz Boosted Z events (MLBZs). The technique is similar to 
the one described in section 4.4.1. The main advantage of this method was that Monte 
Carlo simulated jet properties in W + W~ events could be directly compared with the 
corresponding ones from real data using the large Z statistics. 

The main extension of the method beyond that described in [36] consisted in an 
improved mixing and boosting procedure of the Z events into MLBZs, demonstrated 
in Figure 13. 

The 4-momenta of the four primary quarks in WPHACT generated W + W~ — > qiq~2q3q4 
events were used as event templates. The Z events from data or simulation were chosen 
such that their thrust axis directions were close in polar angle to one of the primary quarks 
of the W + W~ event template. Each template W was then boosted to its rest frame. The 
particles in the final state of a selected Z event were rotated so that the thrust axis matches 
the rest frame direction of the primary quarks in the W + W~ template. After rescaling the 
kinematics of the Z events to match the W boson mass in the generated W + W" template, 
the two Z events were boosted to the lab frame of the W + W~ template. All particles 
having an absolute polar angle with the beam direction smaller than 11° were removed 
from the event. The same generated WPHACT events were used for the construction of both 
the data MLBZs and Monte Carlo MLBZs in order to increase the correlation between 
both emulated samples to about 31%. This correlation was taken into account when 
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Figure 13: Illustration of the mixing and boosting procedure within the MLBZ method 
(see text for details). 
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quoting the statistical uncertainty on the systematic shift on the observables between 
data and Monte Carlo MLBZs. 

It was verified that when introducing a significant mass shift of 300 MeV/c 2 on M w 
by using the cone rejection algorithm, it was reproduced within 15% by applying the 
MLBZ technique. Because the expected systematic uncertainties on the AMw(std,i) 
observables of interest are one order of magnitude smaller than 300 MeV/c 2 , this method 
is clearly justified. 

The double difference of Equation 15 was determined with the MLBZ method using 
Z events selected in the data sets collected during the 1998 calibration runs and Z events 
from the corresponding Monte Carlo samples. The following results were obtained for 
the AM w (std, R C one = 0.5 rad) observable: 



where the statistical uncertainty takes into account the correlation between the Monte 
Carlo and the data MLBZ events, together with the correlation between the two M w 
estimators. This indicates that most of the fragmentation, detector and Between-W 
Bose-Einstein Correlation systematics are small. The study was not performed for the 
AM w (std, p cut ) observable. 

Other systematic sources on the reconstructed jets are not considered as the M w 
estimators used in the difference AM w (std, i) have a large correlation. 

5.2.2 Additional Fragmentation systematic study 

The fragmentation of the primary partons is modelled in the Monte Carlo simulation 
used for the calibration of the Mw* observables. 

The expected values on the M w estimators from simulation (in the k — hypothe- 
sis) are changed when using different fragmentation models [33], resulting in systematic 
uncertainties on the measured Mw* observables and hence possibly also on our esti- 
mated k. In Figure 14 the systematic shift 5Mw in the different Mw* observables is 
shown when using HERWIG or ARIADNE rather than PYTHIA as the fragmentation model 
in the no Colour Reconnection hypothesis. When inferring k from the data difference, 
AM w (std, i), the PYTHIA model is used to calibrate each M w * observable. This data 
difference for M w Pcut=2GeV/c , AM w (std,p cut = 2GeV/c), changes 3 by (27 ± 12) MeV/c 2 
or (8 ± 12) MeV/c 2 when replacing PYTHIA by respectively HERWIG or ARIADNE. Simi- 
larly, the observable AM w (std, R conc = 0.5 rad) changes by (-4 ± 10) MeV/c 2 or (-6 ± 
10) MeV/c 2 when replacing PYTHIA by respectively HERWIG or ARIADNE. The largest 
shift of the observable when changing fragmentation models (or the uncertainty on this 
shift if larger) is taken as systematic uncertainty on the value of the observable. Hence, 
systematic errors of 27 MeV/c 2 for the AMw(std, p cut = 2GeV/c) observable and 10 
MeV/c 2 for the AM w (std, R conc = 0.5 rad) observable were assumed as the contribution 
from fragmentation uncertainties. The MLBZ studies (see above) are compatible with 
these results, hence no additional systematic due to fragmentation was quoted for the 
AMw(std, R C onc = 0.5 rad) observable. 

3 This change, AM w (std, p cut = 2 GcV/c) PYTHIA - AM w (std, p cut = 2 GeV/c) HERWIG , is given by 
SM w (std = p cut = 0.2GeV/c) PYTHIA - HERWIG - 5M w (p cut = 2 GeV/c) PYTHIA - HERWIG , and similar expressions for the 

ARIADNE and Rconc Cases (for Rconc i Std = Rconc = 7r). 



A syst (ARIADNE , DA) 
A sys t (PYTHIA , DA) 
A syst (HERWIG , DA) 



-1.9 ± 3.9(stat)MeV/c 2 
-5.7 ± 3.9(stat)MeV/c 2 
10.6 ± 3.9(stat)MeV/c 2 



(16) 
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Figure 14: Systematic shifts #Mw, on Mw observables, when applying different fragmen- 
tation models as a function of the p cut or R con e values used in the construction of the 
Mw observable. These Monte Carlo estimates were obtained at a centre-of-mass energy 
of 189 GeV. The uncertainties are determined with the Jackknife method. 
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5.2.3 Energy Dependence 

The biases of the different Mw estimators have a different dependence on the centre- 
of-mass energy, hence the calibration of AM w (i, j) will be energy dependent. The energy 
dependence of each individual M w estimator was parameterised with a second order poly- 
nomial. Since WPHACT event samples were used at a range of centre-of-mass energies 
the uncertainty on the parameters describing these curves are small. Therefore a small 
systematic uncertainty of 3 MeV/c 2 was quoted on the AM w (i, j) observables due to the 
calibration. 

5.2.4 Background 

The same event selection criteria were applied for all the M w estimators, hence the 
same background contamination is present in all analyses. The influence of the qq(7) 
background events on the individual Mw estimators is small [33] and was taken into 
account when constructing the centre-of-mass energy dependent calibration curves of 
the individual M w estimators. The residual systematic uncertainty on both AM w (i, j) 
observables is 3 MeV/c 2 . 

5.2.5 Bose-Einstein Correlations 

As for the particle flow method, the systematic uncertainties due to possible Bose- 
Einstein Correlations are estimated via Monte Carlo simulations. The relevant values 
for the systematic uncertainties on the observables are the differences between the ob- 
servables obtained from the Monte Carlo events with Bose-Einstein Correlations inside 
individual W's (BE1) and those with, in addition, Bose-Einstein Correlations between 
identical particles from different W's (BEA). The values were estimated to be (6.4 ± 
9.3) MeV/c 2 for the AM w (std,p cut = 2GeV/c) observable, and (7.2 ± 8.2) MeV/c 2 for 
the AM w (std, R conc = 0.5 rad) observable. As the uncertainties in the estimated contri- 
butions were larger than the contributions themselves, these uncertainties were added in 
quadrature to the systematic errors on the relevant observables. 

5.2.6 Cross-check in the Semi-leptonic Channel 

Colour Reconnection between the decay products originating from different W boson 
decays can only occur in the W + W~ — > qiq2q3q4 channel. The semi-leptonic W + W" 
decay channel (i.e, qq'fX^) is by definition free of those effects. Therefore the determi- 
nation of Colour Reconnection sensitive observables, like AM w (std, R con e = 0.5 rad), in 
this decay channel could indicate the possible presence of residual systematic effects. A 
study of the AM w (std, R conc = 0.5 rad) observable was performed in the semi-leptonic 
decay channel. The semi-leptonic Mw analysis in [33] was used and the cone algorithm 
was implemented in a similar way as for the fully hadronic decay channel. The same 
data sets have been used as presented throughout this paper and the following result was 
obtained: 

AM w (std,R conc ) = M w std - M w Rcono = (8 ± 56(stat))MeV/c 2 (17) 

where the statistical uncertainty was computed taking into account the correlation be- 
tween both measurements. Although the statistical significance of this cross-check is 
small, a good agreement was found for both M w estimators. 



27 



5.3 Results from the Mw Estimators Analyses 

The observable AM w (std, R cone ) with R con e equal to 0.5 rad (defined above), was 
found to be the most sensitive to the SK-I Colour Reconnection model, and the 
AM w (std, p cut = 2 GeV/c) observable was measured as a cross-check. The analyses were 
calibrated with PYTHIA k = WPHACT generated simulation events. The values measured 
from the combined DELPHI data at centre-of-mass energies ranging between 183 and 
208 GeV are: 

AM w (std,R conc ) = M w std - M w Rconc = (59 ± 35(stat) ± 14(syst))MeV/c 2 
AM w (std,p cut ) = M w std - M w Pcut = (143 ± 61(stat) ± 29(syst))MeV/c 2 

(18) 

where the first uncertainty numbers represent the statistical components and the sec- 
ond the combined systematic ones. The full breakdown of the uncertainties on both 
observables can be found in Table 6. 





Uncertainty contribution (MeV/c 2 ) 


Source 


AM w (std, R conc = 0.5 rad) 


AM w (std,p cut = 2GeV/c) 


Fragmentation 


11 


27 


Calibration 


3 


3 


Background 


3 


3 


BEI-BEA 


8 


9 


Total systematic 


14 


29 


Statistical Error 


35 


61 


Total 


38 


67 



Table 6: Breakdown of the total uncertainty on both relevant observables. 



From these values estimates were made for the k parameter by comparing them with 
the Monte Carlo expected values in different hypothesis of k, shown in Figure 15 for the 
observable AM-\y(std, R cone = 0.5 rad). 

The Gaussian uncertainty on the measured observables was used to construct a log- 
likelihood function C(k) = — 21ogL(/t) for k. The log-likelihood function obtained is 
shown in Figure 16 for the first and in Figure 17 for the second observable. 

The result shown in Figure 16 is the primary result of this analysis, because of the 
larger sensitivity of the AMw(std, R cone = 0.5 rad) observable to the value of k (see sec- 
tion 5.1). The value of k most compatible with the data within one standard deviation 
of the measurement is 

«SK-I = l-75l?:S • (19) 
The result on k extracted from the cross-check AMw(std, p cut = 2GeV/c) observable 
is found not to differ significantly from the quoted result obtained with the more opti- 
mal AM w (std, Rcone = 0.5 rad) observable. The significance can be determined by the 
difference between both Mw estimators : 

M w Pcut - M w Rconc = (-84 ± 59(stat))MeV/c 2 . (20) 

Taking into account that the expectation of this difference depends on k, we find a sta- 
tistical deviation of about 1 to 1.5(7 between the measurements. No improved sensitivity 
is obtained by combining the information of both observables. 
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10 _1 1 10 



SK-I Model parameter k 

Figure 15: The dependence of the observable AM\y(std, R CO ne = 0.5 rad) from simulation 
events on the value of the SK-I model parameter k. The dependence is given at three 
centre-of-mass energies. 



29 



o 

CM 



5 
4.5 

4 
3.5 

3 

2.5 
2 

1.5 
1 

0.5 




Likelihood of indirect measurement of k SK-I 

— DELPHI measured : std - R=$ 50 (stat+syst) 

— DELPHI measured : std - R^hf.50 (stat) 




10 



10 10' 
SK-I Model parameter k 



Figure 16: The log-likelihood function —2 logL(/t) obtained from the DELPHI data mea- 
surement of AM w (std, R CO ne = 0.5 rad). The bottom curve (full line) gives the final result 
including the statistical uncertainty on AM w (std, R cone = 0.5 rad) and the investigated 
systematic uncertainty contributions. The top curve (dashed) is centred on the same min- 
imum and reflects the log-likelihood function obtained when only statistical uncertainties 
are taken into account. 
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Figure 17: The log-likelihood function —2 logL(/t) obtained from the DELPHI data mea- 
surement of AM w (std, p cut = 2 GeV/c). The bottom curve (full line) gives the final result 
including the statistical uncertainty on AM w (std, p cut = 2GeV/c) and the investigated 
systematic uncertainty contributions. The top curve (dashed) is centred on the same min- 
imum and reflects the log-likelihood function obtained when only statistical uncertainties 
are taken into account. 
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In this paper the SK-I model for Colour Reconnection implemented in PYTHIA was 
studied because it parameterizes the effect as function of the model parameter k. Other 
phenomenological models implemented in the ARIADNE [7,8] and HERWIG [6] Monte Carlo 
fragmentation schemes exist and are equally plausible. Unfortunately their effect in 
W + W~ — > qiq2q3q~4 events cannot be scaled with a model parameter, analogous to k 
in SK-I, without affecting the fragmentation model parameters. Despite this non- 
factorization property, the consistency of these models with the data can still be ex- 
amined. The Monte Carlo predictions of the observables in the hypothesis with Colour 
Reconnection (calibrated in the hypothesis of no Colour Reconnection) give the following 
values: 



ARIADNE - 


-> M w std 


- M w Rcono 


= (7.2 


± 


4.1) 


ARIADNE - 


-> M w std 


- M w Pcut 


= (9.4 


± 


7.0) 


HERWIG 


-> M w std 


- M w Rcono 


= (19.7 


± 


4.0) 


HERWIG 


-> M w std 


- M w Pcut 


= (22.8 


± 


6.9) 



MeV/c 2 



(21) 



The small effects on the observables with the HERWIG implementation of Colour Reconnec- 
tion compared to those predicted by SK-I are due to the fact that the fraction of events 
that reconnect is smaller in HERWIG (1/9) compared to SK-I (> 25% at y/s = 200 GeV). 
After applying this scale factor between both models, their predicted effect on the W 
mass and on the AMw(ij') observables becomes compatible. The ARIADNE implementa- 
tion of Colour Reconnection has a much smaller influence on the observables compared 
to those predicted with the SK-I and HERWIG Monte Carlo. 



5.4 Correlation with Direct Measurement 

When using a data observable to estimate systematic uncertainties on some measur- 
and inferred from the same data sample, the correlation between the estimator used to 
measure the systematic bias and the estimator of the absolute value of the measurand 
should be taken into account. Therefore the correlation between the Colour Reconnection 
sensitive observables AM w (std, R con c = 0.5 rad) and AM w (std, p cut = 2GeV/c) and the 
absolute M\y(std) estimator was calculated. The correlation was determined from the 
Monte Carlo events and with k = or no Colour Reconnection. The values obtained 
were found to be stable as a function of k within the statistical precision. The correlation 
between AM w (std, R con e = 0.5 rad) and M w (std) was found to be 11%, while for the one 
between AM w (std,p cut = 2GeV/c) and M w (std) a value of 8% was obtained. Also the 
correlation between the different Mw estimators was estimated and found to be stable 
with the value of k. A value of 83% was obtained for the correlation between Mw(std) 
and M w Rconc=a5rad , while 66% was obtained between M w (std) and M w Pcut=2GcV/c . 

6 Combination of the Results in the Scope of the 
SK-I Model 

The log-likelihood curve from the particle flow method was combined with the curve 
from the AM W method and the result is shown in Figure 18. The correlations between the 
analyses were neglected because the overlap between the samples is small and the nature 
of the analyses is very different. The total errors were used (statistical and systematic 
added in quadrature) in the combination. 
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Figure 18: The log-likelihood function —2 logL(/t) obtained from the combined DELPHI 
measurement via AMw(std, R conc = 0.5rad) and the particle flow. The full line gives 
the final result including the statistical and systematic uncertainties. The log-likelihood 
functions are combined in the hypothesis of no correlation between the statistical and 
systematic uncertainties of both measurements. 
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The best value for k from the minimum of the curve, with its error given by the width 
of the curve at the value — 21ogL = (— 21ogL) min + 1, is: 



«SK-I = 2.21?:! • (22) 



7 Conclusions 

Colour Reconnection (CR) effects in the fully hadronic decays of W pairs, produced 
in the DELPHI experiment at LEP, were investigated using the methods of the particle 
flow and the Mw estimators, notably the AMw(std, R con c = 0.5 rad) observable. 

The average of the ratios R of the integrals between 0.2 and 0.8 of the particle distri- 
bution in Inside-W regions to the Between-W regions was found to be 

(R) = 0.979 ± 0.032(stat) ± 0.035(syst) . (23) 

The values used in this average were obtained after rescaling the value at each energy to 
the value at a centre-of-mass energy of 189 GeV using a fit to the MC without CR. 

The effects of CR on the values of the reconstructed mass of the W boson, as imple- 
mented in different Monte Carlo models, were studied with different estimators. From 
the estimator of the W mass with the strongest sensitivity to the SK-I model of CR, the 
AM w (std, Rconc = 0.5 rad) method, the difference in data was found to be 

AM w (std, R cone ) = M w std - M w Rconc=0 - 5rad = ( 59 ± 35(stat) ± 14(syst) )MeV/c 2 . (24) 

From the combination of the results from particle flow and M-yv estimators, corre- 
sponding to the curve in full line shown in Figure 18, the best value and total error for 
the k parameter in the SK-I model was extracted to be: 

K SK-I = 2-2±?:| (25) 

which corresponds to a probability of reconnection of V rcC o = 52% and lies in the range 
31% < V TCCO < 68% at 68% confidence level. 

The two analysis methods used in this paper are complementary: the method of parti- 
cle flow provides a model-independent measurement but has significantly less sensitivity 
towards the SK-I model of CR than the method of AM W estimators. 

The obtained value of k in equation (25) can be compared with similar values obtained 
by other LEP experiments, and it was found to be compatible with, but higher than, the 
values obtained with the particle flow by L3 [37] and OPAL [38]. It is also compatible 
with, but higher than, the values obtained with the method of different Mw estimators 
by OPAL [39] and ALEPH [40]. 
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